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ABSTRACT 


Two algorithms are used to determine bathymetry in the littoral region using 
QuickBird multi-speetral satellite observations. The algorithms determine water-leaving 
radianee and eonvert this to water depth values. The first algorithm uses a ratio of two 
wavebands and the seeond uses the sum of several wavebands. Relative bathymetrie 
errors are determined for the elear water of Looe Key (USA) and the turbid water of 
Plymouth Sound (UK). Bathymetrie measurements from LIDAR and ehart data are 
eompared to derived depths to assess their aceuracies. An amended version of the ratio 
method is proposed for use in turbid water to improve aeeuraey. The results show that 
the standard ratio and turbidity algorithms have a relative error of 11.7% and 16.5% 
respeetively in elear water. In turbid water the average error of the turbidity algorithm is 
11.6% and the amended ratio algorithm average error is 13%. 
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I. 


INTRODUCTION 


Remote sensing has beeome an inereasingly important applieation for bathymetrie 
mapping sinee its ineeption as a tool to determine water depth over twenty years ago. As 
a result of an inerease in the availability of commereially provided high-resolution multi- 
speetral satellite imagery, such as Quick Bird, SPOT, IKONOS and Landsat, the accuracy 
of water depth analysis has improved accordingly (Stumpf et al, 2003). Given that 
bathymetry can be determined at substantially less cost and more expediently by remote 
sensing techniques than traditional hydrographic methods, an analysis of its accuracy in a 
range of environments is proposed to determine usability in differing operational 
environments. 

Remote sensing data can potentially service a range of requirements including 
navigation data, updates to the Digital Bathymetric Data Base (DBDB), inputs to 
modeling algorithms for surf, current and acoustic predictions, mine detection. Beach 
Intelligence Survey Data and environmental monitoring. In order to achieve these 
requirements, a degree of accuracy must be achieved. Current accuracy comparisons 
with hydrographic data show mean errors in the range of 10-30% with a maximum 
effective depth ranging from 6-10.5m (Lafon et al, 2002, Tanis et al, 2002). The major 
source of error is the compounded attenuation by the atmosphere and water column on 
radiance reflected by the seabed and received by the satellite. 

The basic premise behind water depth calculations is that the radiance received by 
the satellite is a function of five parameters: incoming solar radiation; attenuation of 
radiation into and out of the atmosphere; attenuation of radiation into and out of the water 
column; reflectance properties of the seabed and the depth of water. The depth of water 
can be resolved by finding values for the first four parameters (Lyzenga, 1981). Several 
techniques exist to achieve this, they fall into three main categories: 

1) The first are linear methods that solve empirically for several variables, 
then use Beer’s Law, which states that light decays exponentially with 
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depth in the water column (Lyzenga, 1981, 1985). A single waveband is 
used to determine water depth for each pixel. 

2) The second are ratio methods, which use two or more wavebands to 
determine the depth of water based on the differing attenuation properties 
of the multiple wavelengths. Radiance values from these wavebands solve 
for depth through correlation with known ratios for a specific water 
column. These methods have fewer parameters to solve for and are less 
affected by changes in bottom reflectance as both wavelengths are affected 
similarly by bottom albedo (Stumpf et al, 2003). 

3) The third type consists of variations on the linear method. The Stratified 
Genetic Algorithm uses a method that splits the water column into levels 
of differing attenuations (due to turbidity) and calculates a sum of ratios 
for several wavelengths (Gianinetto et al, 2003) 

The potential for error exists because of the difficulty in accurately determining 
the optical properties of the atmosphere, ocean and seabed. Unfortunately model data 
must be used as there is little chance in an operational situation to take in situ 
measurements at the time of image capture. The objectives of this thesis are to compare 
each method’s performance in a range of different environmental conditions and evaluate 
relative accuracies against hydrographic and Lidar survey data. 

Chapter II provides a review of relevant atmosphere and ocean optical properties 
and the radiative transfer equation. Chapter III outlines the data used. Chapter IV 
discusses the applications and methods used to analyse the data. Chapter V presents the 
results and Chapter VI summarises the conclusions and recommendations. 
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II. 


BACKGROUND 


Central to retrieving water depth information from remotely sensed radianee 
values are aeeurate eorreetions for attenuation and baek-scattering within the atmosphere 
and water column. Radiative transfer theory relates the measured radiance at the sensor 
to the interaction between incoming solar radiation, atmospheric aerosols (Martin, 2004) 
and water column attenuation (Jupp, 1988). Corrections for the effects of the atmosphere 
and water column must be applied to the remotely sensed data to isolate the measured 
radiance as a function of water depth. The total radiation measured at the satellite is the 
product of incoming solar radiation, attenuation within the atmosphere and water column, 
atmospheric backscatter, water surface reflectance, substrata reflectance and water depth. 



Eo Solar radiance at the top of 

atmosphere 

Ed Down-welling irradiance 

Ea Diffuse sky irradiance 

Lws Reflected light from water 

surface 

Rw Water column irradiance 

reflectance 

Ad Bottom irradiance 

reflectance 

R(0-) Irradiance reflectance just 

below water surface 
Lw Water-leaving radiance 

Eu Upwelling irradiance 

Lp Atmospheric path radiance 

Lp Radiance transmitted from 

the top of the atmosphere 
Ed Radiance received at the 

sensor 

Ta Atmospheric transmittance 


Figure 1. Factors that influence radiance reaching a sensor over a water mass. (After: 

Bierwirth et al,1992). 
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A. OPTICAL PROPERTIES OF THE ATMOSPHERE 

Radiative transfer theory outlines four main sourees of photons that eontribute to 
the total radianee reeeived by a sensor at the top of the atmosphere (Durkee et al, 1986); 

L,=L, + L.+(l,+L.)t„ ( 2 . 1 ) 


where Lt is the total radianee reeeived by the sensor, Lr is the moleeular Rayleigh 
seattering radianee, La is the aerosol seattering radianee, Lg is the radianee from surfaee 
glint, Ls is the surfaee radianee and Td is the transmittanee of the of the direet path. 
Atmospherie seattering is ealeulated to separate the surfaee radianee from the total 
radianee. A surfaee reflectanee algorithm is used to ealeulate the amount of refleetion 
produeed by surfaee glint (or speeular reflection) and this value is then removed from the 
top of the atmosphere radiance. Model calculations are used to assess molecular 
Rayleigh scattering within different wavebands and the information incorporated into the 
NPS aerosol model. Using linearized single scatter theory (LSS), this model solves for 
overall atmospheric scattering radiance given as: 






( 2 . 2 ) 


where k is wavelength, 0 is the satellite zenith angle, (p is the satellite azimuth angle, Lq is 
surface radiance, is the extinction term, ©o is the single scatter albedo, p(v|/s) is the 

scattering phase function, 5 is optical depth, 0o is the sun zenith angle, Oq is the sun 
azimuth angle and p is cosine of the satellite zenith angle (0). Solving for the radiance 
due to aerosol scattering becomes: 


L„ = 


a. 




P{¥s)Sa^ 


(2.3) 


where 6a is aerosol optical depth. This reduced transfer equation is based on the 
following assumptions: the region is sun glint free, the aerosols are spherical and non¬ 
absorbing and single scattering dominates (Durkee et al, 1991). This model was applied 
to the QuickBird imagery to determine aerosol optical depth (AOD). 
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B. OPTICAL PROPERTIES OF THE OCEAN 

Knowledge of the radianee within and leaving the water eolumn is a prerequisite 
to deriving a depth value from Lt eorreeted for atmospherie attenuation and seattering. If 
the optieal properties of water are known and bottom refleetanoe is uniform, a good 
eorrelation between water depth and a waveband ean be assumed (Lyzenga, 1985). 

Several eonstituents of water signifieantly alter the seattering and attenuation of 
radiation in the visible and near-infrared speetrum. Salts inerease seattering by 30% in 
eomparison to fresh water but have little effeet on absorption. Raman seattering 
generates many different wavelengths of light from a nominally single-wavelength souree 
through interaetion with moleeules within the water eolumn. Organie eompounds (or 
yellow matter) and Coloured Dissolved Organie Matter (CDOM) are ereated by deeaying 
plant matter. CDOM is a signifieant absorber of the blue waveband. Inorganie 
partieulate matter sueh as fine minerals from estuarine outlets ean seatter visible light. 
Organie particulate matter, including chlorophyll, viruses, bacteria, phytoplankton and 
organic detritus can, depending on their concentration and size, scatter visible light 
appreciably (Fisher, 1999). 

The surface morphology of the sea surface also affects the radiance of light from 
the water column. Ripples and waves alter the reflective properties of the sea surface and 
introduce variability, or noise, in the received radiance values at the sensor (Clark, 2005). 


5 




500 600 700 800 900 


Wavlength (nm) 

Figure 2. Water refleetanee speetra measured by speetro-radiometer with varying chlorophyll 

amounts. (From: Lavender et al, 2004). 


The sum of spectral absorption and scattering coefficients are termed the Inherent 
Optical Properties (lOPs) of water and are related to the constituents of the water column. 
The spectral absorption coefficient is the fraction of incident power at a given wavelength 
that is absorbed per unit distance in a medium. The spectral scattering coefficient is the 
amount of incident power per unit distance that is scattered out of the beam (Mobley, 
1995). Whereas lOPs are dependent on the medium alone, Apparent Optical Properties 
(AOPs) are dependant on the directionality of radiance. The spectral radiance reflectance 
(Rs) is the ratio of spectral upwelling irradiance to spectral downwelling irradiance at a 
given depth: 




EAz,X) 


(2.4) 


where z is depth, X is wavelength, Eu is upwelling irradiance and Ed is downwelling 
irradiance. Spectral remote sensing reflectance (R) is the ratio of water leaving radiance 
(Ew) to downwelling irradiance (Ed) and is calculated just below the surface: 


£„(o-^) 






E,(0-,^)’ 


(2.5) 
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where (0-) denotes some small depth just below the water’s surface. This value is the 
amount of downwelling light that leaves the surface after interacting with the water 
column and the seabed and is recorded by the remote sensor (Fisher, 1999). 

C. THE RADIATIVE TRANSFER EQUATION 

To determine water depth from satellite-recorded radiances, atmospheric and 
water column radiative transfer properties are summed and the result is solved for water 
depth. There are three main methods: 

1. Linear Method 

Linear methods use the assumption that light attenuates as a power law in water to 
determine depth (Lyzenga, 1978; Jupp, 1988; Philpot, 1989) and that this attenuation is 
spatially homogeneous. Philpot (1989) states: 

L, = , (2.6) 

where La is the radiance observed at the remote sensor. Lb is the bottom reflected 
radiance, g is the water attenuation coefficient, z is depth and Lw is the radiance over 
optical deep water. The assumption is that atmospheric transmittance, path radiance, sea 
state, water reflectance, bottom reflectance, and water attenuation are constant. In order 
to incorporate atmospheric and water column optical properties, this formula is 
rearranged to give: 

L, (z) = CE, {0-){Ab - p + CE, (0-)p„ + , (2.7) 

where C is the atmospheric transmission factor, Ed(O-) is downwelling irradiance just 
below the water surface, Ab is bottom reflectance (albedo), poo is the irradiance 
reflectance of optically deep water, Ljg is the sun glint and Lpath is the path radiance. The 
atmospheric and sun glint terms are removed and the radiance values are converted to 
top-of-the-atmosphere values. Melsheimer and Liew (2001) state: 

R(z) = ^R(0)-R(oo)]-^^+R(ao), (2.8) 

where R is top-of-the atmosphere reflectance. Equation 2.8 is then solved for depth: 
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z = - (log[i?(z) - i?(G0)] - log[i?(0) - i?(G0)]), 


(2.9) 


2. Ratio Method 

Ratio methods build on the assumption that light attenuates exponentially with 
water depth by comparing the ratio of different attenuation rates between wavebands to 
determine depth. As depth increases, the water leaving radiance of a band with higher 
absorption (green) will decrease proportionately faster than that of a band with lower 
absorption (blue) and the ratio of blue to green values will therefore increase (Stumpf and 
Holderied, 2003). The basis of the ratio transfer equation is the linear relationship 
presented by Lyzenga (1985) that states: 

Z=a,X.+ajXj+Z^, ( 2 . 10 ) 

with: 

X, =ln[/;a,)-S„(.l,)], (2.11) 


where R is the reflectance at wavelength (>.i), Roo is the reflectance of optically deep water 
at wavelength (ki), and ai, aj and Zq are coefficients that account for the optical properties 
of the water. This is combined as a ratio of two wavelengths by Stumpf and Holderied 
(2003) to give: 


Z = mj 


ln{nRJAj)) 


( 2 . 12 ) 


where mi is a tunable constant to scale the ratio to depth, n is a fixed constant for all areas 
to ensure the terms remain positive and mo is an offset for a reference depth of Om (the Zo 
term in equation 2.10). The theoretical benefit of a ratio transform is a compensation for 
variable bottom reflectance. It is claimed a change in bottom albedo affects both bands 
similarly, but a change in depth affects the higher waveband more, so variable bottom 
types can be analysed more accurately. 


3. Stratified Genetic Algorithm 

The Stratified Genetic Algorithm (SGA) is a development of the Depth of 
Precision (DOP) model proposed by Jupp (1988) that states: 
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(2.13) 


Le is measured radianee at the sensor, Lb is radianee from the seabed, k is the absorption 
eoeffieient of the water and N is the number of speetral bands. The seeond term is 
removed and replaeed with a regression eoeffieient (Yj) by Gianinetto and Leehi (2003) 
to give: 


f.ln(4).-T. 


(2.14) 


where m is the number of layers. The SGA method divides the water eolumn into levels 
of inereasing depth and eomputes kj and Yj for eaeh in order to ealeulate water depth. 
This algorithm is repeated for all speetral wavebands and those with a high eorrelation 
eoeffieient are used to determine depth. 
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III. DATA 


A. DATA DESCRIPTION 

The data used ineludes Digital Globe’s QuiekBird (QB) polar orbiting satellite 
image retrieval system, LIDAR bathymetry soundings, United Kingdom Hydrographie 
Offlee (UKHO) and National Geophysieal Data Centre (NGDC) hydrographie ehart data, 
and the HYDROLIGHT water eoeffieient model. 


B. INSTRUMENTS 

I. QuiekBird Data 

The QB 2 satellite is a polar-orbiting, sun-synehronous satellite operating at an 
altitude of 450 kilometers with a 98-degree inelination on a 93.4 minute orbit eyele with a 
10:30 a.m. loeal equator erossing time for the aseending node (Digital Globe, 2004). The 
satellite eaptures solar refleeted energy in four speetral bands in a swath 16.5km wide. 


CHANNEL 

BANDWIDTH (nm) 

SPATIAL RESOLUTION(m) 

1 

450-520 

2.44-2.8 

2 

520-600 

2.44-2.8 

3 

630-690 

2.44-2.8 

4 

760-900 

2.44-2.8 


Table 1. QB bandwidth and maximum spatial resolutions (DigitalGlobe 2004). 

The multispeetral bands shown in Table 1 have a nominal spatial resolution of 
2.44-2.8m at nadir and in addition a panehromatie band (not shown) that has a nominal 
resolution of 0.6m meters. The QB data was aequired by the National Geospatial Ageney 
(NGA) after being proeessed and geo-refereneed by DigitalGlobe ground stations. 
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Catalog ID: 101001000382D601 Acq Date; Jan. 13 2005 LaVLong: 50.3324* / •4.16485* 

Off NadirAngle: 9* TargetAzimuth; 350* Cloud Cover: 3% Environmental Quality; 70 


DiGITALGLOBE 

<e>r)iait«)iaiob» zoos ah Riaht» r 


To Order Call: (303}6S4~4561 or (800)49S-1225 
Email: orders@digitalglot>e.com 
Or Click Here To Find Your DigitalGlobe Reseller 


Figures. Digital Globe’s QB imagery of Plymouth Sound. (After; DigitalGlobe, 2005). 

2. LIDAR Data 

Light Detection and Ranging (LIDAR) determines range records using the time 
taken for an emitted pulse of energy to hit a target and return to the sensor to determine 
range. The Scanning Hydrographic Operational Airborne LIDAR Survey (SHOALS) 
system uses this technology to determining water depth. SHOALS is managed by the 
Joint Airborne LIDAR Bathymetry Technical Center of Expertise (JALBTCX). The 
airborne component of the system is installed aboard a DHC/300 Twin Otter airframe and 
uses a 1064nm primary pulse for an infrared spectra and a 532nm pulse for a blue-green 
spectra to collect bathymetry data on a 4m grid. 
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MAXIMUM DEPTH 

60 m or 2-3 times Secchi depth 

VERTICAL ACCURACY 

+/- 15 cm 

HORIZONTAL ACCURACY DGPS 

+/- 2 m 

HORIZONTAL ACCURACY OTF GPS 

+/- 0.5 m 


Table 2. SHAOLS performance specification (After; http;//shoals.sam.usace.army.mil, 

used 15 Mar 2005). 

LIDAR soundings were acquired for Looe Key on 10 December 2004 by 
JALBTCX and projected to the North American 1988 Datum. The depth of the seabed 
was calculated below the water surface at the time of data capture. A height of tide 
correction and datum conversion to WGS 84 is made before comparison with satellite 
derived depths. 

3. UKHO and NGDC Hydrographic Data 

Satellite depth analysis is also compared to navigation chart data supplied by the 
UKHO for the Plymouth and the NGDC for the Looe Key area. These products comply 
with the International Hydrographic Organisation (IHO) standards for hydrographic 
surveys of Order 1. This order is required for harbours, harbour approach channels, 
recommended tracks, inland navigation channels, and coastal areas of high commercial 
traffic density (less than 100m) with depths up to 100 m (S44 IHO 1998). The basic 
assumption is that errors are classified by: a, the sum of all dept-independent errors and b, 
the sum of all depth-dependent errors, expressed as a percentage of water depth. The 
degrees of accuracy required are divided into four categories depending on the product to 
be manufactured. The inshore navigation charts used in this study comply with the 
Category 1 specification outlined in Table 3. 
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ORDER 

SPECIAL 

1 

2 

3 

Depth of uncertainty for 

a = 0.25m 

a = 0.5m 

a = 1.0m 

a = 1.0m 

reduced depths (95% 
confidence level) 

b = 0.75% 

b= 1.3% 

b = 2.3% 

b = 2.3% 

Bathymetric model 

a = 0.25m 

a = Im 

a = 2m 

a = 5m 

uncertainty (95% 
confidence level) 

b = 0.75% 

b = 2.6% 

b = 5% 

b = 5% 


Table 3. Summary of Minimum Standards for Depth Uncertainties. (After; S44 IHO 

1998). 


The UKHO charts are projected to the WGS 84 datum, whereas the NGDC chart 
is projected onto the North America Datum of 1983. NGDC depths were corrected to 
WGS 84 before depth comparisons with QB derived depths were made. Because QB 
imagery is processed to WGS 84, the NGDC depths were corrected to WGS 84 before 
depth comparisons with satellite derived depths were made. 



Figure 4. Nautical charts for Plymouth Sound and Looe Key (After; UKHO, NGDC). 
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4. HYDROLIGHT 

The HYDROLIGHT radiative transfer numerieal model eomputes radianee 
distributions and related quantities ineluding irradianees, refleetanees and diffuse 
attenuation funetions of the water eolumn. Inputs from the user inelude the water 
absorption and seattering properties, the sun and sky radianees, and the bottom 
eharaeteristies. The output ineludes the water-leaving radianee and remote-sensing 
refleetanee of the water eolumn, known as ^-funetions. To solve for water depth, the 
downwelling attenuation eoeffieient output {Kd), is determined for the QB wavebands 
channels used and incorporated into the various depth algorithms. 

A classification scheme has been developed for oceanic waters based on the value 
of Kd (Jerlov, 1976). The three different optical types. Case I, Case II, and Case III 
classify the increase in attenuation, scattering and absorption of photons with differing 
levels of turbidity. Case I (divided into lA and IB) are the clearest and Case III being the 
most turbid. The HYDROLIGHT 3.0 version is used to determine Kd for Case I and Case 
II waters. The Case II waters are used for Plymouth Sound with an observed value for 
chlorophyll concentration, via a separate FORTRAN sub-routine, and supplied values for 
CDOM and Raman scattering. Case I water is used for Looe Key where CDOM and 
particulate matter within the water column is greatly reduced. An averaged value for 
chlorophyll is used for each area of interest as the lengthy run times for evaluating Kd 
prohibits evaluating the whole image for a unique value for each pixel. 
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IV. METHODOLOGY 


A. ATMOSPHERIC CORRECTION 

The amount of energy reeeived by a remote sensor from atmospherie aerosol 
effeets and surfaee refleetanee adds to the amount of energy radiated from the water 
eolumn and seabed and must therefore be ealeulated and eorreeted for. Aerosol 
eharaeterisation information is aequired from in situ or remote methods. In situ 
measurements are made using ground or air based radiometers, sun-photometers and 
speetrometers, and provide high resolution data for highly defined areas. Remotely 
sensed methods, sueh as those used with the Sea-viewing Wide Field-of-view Sensor 
(SeaWIFs) remote sensing system, use satellite mounted radiometers that provide an 
analysis over a large spatial area at lower resolution. Without aeeess to in situ 
measurements, aerosol optieal depth (AOD) of the atmosphere is ealeulated from the QB 
imagery using the over-water dark objeet approaeh and the NFS Aerosol Model. Sun 
glint and other refleetanee effeets are eorreeted for using a near infrared (NIR) sealing 
algorithm applied to the QB imagery. 

I. NPS Aerosol Model 

The NPS Aerosol Model applies linearised single seatter theory to an estimate of 
bidireetional surfaee refleetanee using seattering phase funetions that are parameterized 
into seven aerosol model size distributions (Durkee et al, 1991). The aerosol model sizes 
are based on those typieally found in a marine environment and a partiele size parameter, 
ealled S12, is ealeulated using the ratio of the red and NIR wavelength radianees. The 
S12 value and seattering angle are used to determine aerosol model index (AMI), and 
Mie theory is used to ealeulate seattering phase funetion curves for each aerosol model 
size distribution. The AMI and scattering angle determine the scattering phase function 
that is needed to calculate AOD. 

The NPS Aerosol Model was designed for use with NOAA-14 AVHRR radiance 
values. In order to apply the model to the different wavebands of the QB data, the 
radiance values are adjusted. The NOAA-14 AVHRR radiances were linearly 
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extrapolated to rnateh the effeetive eenter red and NIR wavelengths of the QB 
wavebands. The linear extrapolation method is eonsidered parsimonious due to the linear 
nature of in-band solar irradianee and nearly linear nature of the seattering meehanisms at 
the two wavelengths. The extrapolated model radianees were used in the NFS Aerosol 
Model as the input for top of the atmosphere radianee. The output values for QB red and 
NIR eentral wavelengths for the respeetive wavebands were then extrapolated linearly to 
the blue, green and red waveband eentral wavelengths. The retrieved aerosol optieal 
depths for eaeh area are presented in Table 4. 


LOCATION 

BLUE AOD 

GREEN AOD 

RED AOD 

NIR AOD 

Flymouth Sound 

0.3778 

0.2311 

0.1468 

0.0090 

Looe Key 

0.3019 

0.2617 

0.231 

0.1028 


Table 4. NFS Aerosol Model values of AOD for QB imagery. 


2. Sea Surface Correction 

The high resolution of QB imagery inereases the ineidenee of quasi-stoehastie sea 
surfaee effeets, eaused by ripples and waves, whieh ereates speeular refleetanee. This aets 
as a signifieant souree of additional energy reeeived at the remote sensor. An inereased 
surfaee wind results in greater seattering and ereates glint that partially obseures the 
water leaving radianee that is fundamental in determining water depth. The Sea Surfaee 
Correetion algorithm (Hoehberg et ah, 2003) is used to eliminate the majority of the wave 
and glint effeets. This model utilises the NIR waveband, whieh exhibits the maximum 
absorption and minimal water leaving radianee over elear waters, to eharaeterise the 
spatial distribution of relative glint intensities. It is assumed that although absolute glint 
intensity varies with wavelength, relative glint intensity is eonstant aeross NIR and 
visible wavelengths. The NIR eorreeted glint distribution is sealed in the blue, green and 
red QB wavebands, and the respeetive results are subtraeted from the measured radianees 
to remove glint effeets. Uneorreeted and eorreeted images for Looe Key ean be seen in 
Figs. 5 and 6. 
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Figure 5. Image of Looe Key before sea surfaee eorreetion was applied showing ripples on the 

sea surfaee. 





3. Conversion to Reflectance Values 

After removing the aerosol and glint effeets, the resultant QB imagery is 
converted from radiance values to reflectance values using Eq. 4.1 for use in the various 
depth retrieval algorithms; 


Eo 

where R is reflectance, Rs is surface radiance value and Eo is the incoming solar 
irradiance. 

B. WATER COLUMN CORRECTION 

The effective attenuation coefficient of the water mass must be determined in 
order to relate the water depth to the reflectance values received at the sensor. The 
HYDROEIGHT model is used to determine a value for the attenuation coefficient based 
on a value of chlorophyll determined using the QB imagery. 

1. Chlorophyll Analysis 

The amount of chlorophyll in each image is calculated by determining the relative 
chlorophyll levels and converting this to a concentration figure. Radiance at all 
wavelengths decreases with an increase in chlorophyll concentration and this rate of 
decrease differs between wavebands. The rate of decrease is more rapid for blue than 
green wavelengths. A relative value for chlorophyll concentration can be derived from a 
ratio of the reflectance at the blue and green wavelengths. The spatial distribution of the 
chlorophyll concentration can be seen for each area in Eigs. 7 and 8. 
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Figure 7. Looe Key ratio of blue/green wavebands (relative ehlorophyll eontent) 

derived from QB imagery. 



Figure 8. Plymouth Sound ratio of blue/green wavebands (relative ehlorophyll eontent) derived 

from QB imagery. 
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The ratio values are then be eonverted to a eoneentration amount using a pigment 
algorithm for Case 1 and 2 waters (Gordon and Clark, 1981): 

Chl=A^ (4.2) 

where Chi is the ehlorophyll eoneentration (mg/m ), A and B are eonstants, ki and kj are 
the blue waveband radianee and green waveband radianees respeetively. The variability 
in eoneentration values aeross eaeh study area is resolved by seleeting one area in eaeh 
image where depth ealeulations are to be earried out and an average value determined for 
input into the HYDROLIGHT model. The retrieved ehlorophyll eoneentrations are 
shown in Table 5 for the areas shown in Figs 7 and 8. 


AREA 

CHL (mg/m^) 

Plymouth Sound 

0.303 

Looe Key 

0.266 


Table 5. Caleulated ehlorophyll eoneentrations for study areas. 

2. HYDROLIGHT 

The wavelength dependent attenuation eoeffieient, IQ, is required to derive depth 
using the linear and the stratified genetie algorithms. These were modeled using 
HYDROLIGHT 3.0 and ineorporating pre-delined settings for Case 1 water for Looe 
Key and Case II water for Plymouth Sound. A separate FORTRAN subroutine was used 
to fix the value of ehlorophyll at the given eoneentrations for eaeh area, as 
HYDROLIGHT would allow the ehlorophyll eontent to vary with depth otherwise. 

The HYDROLIGHT input parameters were average ehlorophyll eoneentration, 
pressure, temperature, humidity, visibility, preeipitation, longitude, latitude, wavelength, 
solar zenith angle and satellite zenith angle. The output presented in Table 3 shows that 
as wavelength inereases so does Kd, as would be expeeted due to higher absorption of 
visible light by water at longer wavelengths. The ehlorophyll and CDOM amounts at 
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Plymouth Sound were significantly higher than Looe Key, and this is reflected in the 
higher attenuation values at these locations. 

The HYDROLIGHT results were compared with empirically observed Kj values. 
These were derived from re-arranging the linear depth algorithm at known depths to 
solve for IQ. The HYDROLIGHT results are shown in Fig 9. 


Looe Key Water Attenuation 



Figure 9. HYDRO LIGHT derived water attenuation coefficient (IQ) for each area. 


C. DEPTH RETRIEVAL 

An analysis of the stratified genetic and linear algorithm revealed that both 
models were based on the same derived form of the radiative transfer equation. The 
stratified genetic algorithm differs from the linear algorithm by incorporating a regression 
analysis-derived coefficient for the second term of the linear algorithm; 
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z = — (log[i?(z) - R{cc)] - log[i?(0) - i?(oo)]) 


(2.9) 


z = 


P - 2 ^. 


(2.14) 


where kj the water attenuation coeffieient and g is 2kj, Le and R(z) are the measured 
refleetance values and Yj, and R(oo) are the measured reflectance values over optically 
deep water. The coefficient Yj was determined by empirical analysis using known depths 
along each transect. The stratified genetic algorithm was used in comparison with the 
ratio algorithm to retrieve depths for each area. 


1. Ratio Transfer Algorithm 

The ratio algorithm requires empirically derived coefficients to normalise the 
natural log of the ratio of green to blue reflectance values to retrieve a depth term and 
then reference those depths to the sea surface. The first coefficient, mi_ is a tunable 
constant to scale the ratio to a relative depth in meters. The second coefficient, m 2 , is an 
offset value to reference the first term to Om depth. The ratio algorithm assumes that “a 
change in bottom albedo affects both bands similarly, but changes in depth affect the high 
absorption band more”.i A regression algorithm was used with two known depths from 
each area to determine respective values of mi and m 2 respectively. 

In the turbid Plymouth Sound waters, the blue waveband was scattered to such an 
extent that valid radiance values were not possible for use in the ratio algorithm. To 
resolve this problem, the radiance levels for each waveband were analysed along a 440m 
transect in the estuarine waters of the Sound. The results, detailed in Figs 10-12, show 
that the blue waveband is scattered significantly and no discemable correlation can be 
made between an increase in depth and reduction in radiance value. The green waveband 
shows a degree of decay in value with increased depth along the transect, and the red 
band shows the greatest amount of response to depth. The effect of scattering by CDOM 


1 Determination of water depth with high-resolution satellite imagery over variable bottom types. 
Stumpf and Holderied. Ameriean Society of Limnology and Oceanography, 2003. 


24 



and particulate matter in the red and green waveband is less than in the blue, so therefore, 
it is proposed that an amendment to the ratio technique is employed in turbid waters, 
namely, the use of the green to red reflectance values. In this analysis the standard ratio 
method was used for Looe Key and the amended version for Plymouth Sound. 


Blue Radiance Values (Plymouth Transect) 



Figure 10. Digital Number (DN) radiance values for QB blue waveband. 
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2. Stratified Genetic (Turbidity) Algorithm 

The second term of the genetic algorithm, Yj, was derived through a regression 
analysis in a similar manner to the ratio method. Bottom reflectance values for each 
wavelength were taken from the QB images where the wet bottom type was exposed, or 
lay a few centimeters beneath the water’s surface. The area for analysis in each image 
featured a predominantly sand bottom type. This value was incorporated into a term that 
included bottom albedo and the reflectance value of optically deep water: 




-RA 

n _ 


(4.3) 


where Rb is bottom reflectance and Rd is reflectance of optically deep water. 


3. Environment for Visualising Images (ENVI) 

The Research Systems Incorporated (RSI) ENVI software package was used to 
orthorectify and calibrate the QB imagery for accurate comparison with the LIDAR and 
chart data. The Band Math tool was used to execute the correction algorithms for aerosol 
optical depth and sea surface reflection and to convert radiance values to reflectance 
values for each pixel. The two algorithms were then entered into Band Math and 
executed for each image to calculate depth values for each pixel. Transects were taken 
within each area for cross comparison with LIDAR and hydrographic chart data. 
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V. RESULTS 


A. TEST AREAS 

To compare the results over a variety of bottom types and depths, a transect was 
selected within each area, and both the ratio and turbidity algorithms are used to 
determine water depth along the transect. Accuracy was calculated as the root mean 
square (rms) error of the derived depths from LIDAR and chart data. 


1. Looe Key 

The transect at Looe Key extends for 646m in a south-easterly direction from the 
shallow waters of the key to the edge of a submarine drop off. The depths range from 
3.5m at the shallower northwestern end of the transect, increasing to 14.17m towards the 
drop-off. The bottom types encompass sea grass at the northern end of the line with dark 
coral fingers roughly two thirds down. LIDAR data is used to assess accuracy with the 
remotely sensed derived depths. The test transect extending from 24“ 32’ 57.53”N, 81“ 
24’ 10.55”W to 24 32 4L17N, 81“ 23’ 55.99”W is shown in Fig 13: 




2. Plymouth Sound 

The transect in Plymouth Sound extends 1115.3m in a south-southwesterly 
direction from the southern shore of Mount Batten Point into the Tamar estuary. The 
depths range from 0.5m inshore to 10.6m in the estuary. The bottom type was assumed 
to be sand throughout the area of interest. A third of the way along the transect, a visible 
band of more sedimentary water was identified. The transect was split into a smaller 
368m length that marked the extent of the cleaner inshore water with which accuracy 
comparisons with chart data were made. The effect of the more sedimentary water was 
also compared to charted depths. The long test transect extending from 50° 21’ 30.12”N, 
004° T 47.01”W to 50° 20’ 54.28”, 004° 7’ 53.64”W is shown in Fig 14: 



Figure 14. Plymouth Sound test transect. 
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B. DERIVED DEPTHS 


1. Looe Key 

The QB derived depths using the ratio and turbidity algorithms in comparison 
with the LIDAR bathymetry are shown in Figs 15 and 16; 


Looe Key Calculated Depths (Ratio Algorithm) 



Figure 15. Ratio algorithm and LIDAR bathymetry of Looe Key transect. 



Looe Key Calculated Depths (Turbidity Algorithm) 
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Figure 16. Turbidity algorithm and LIDAR bathymetry of Looe Key transect. 
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2. Plymouth Sound (Short Transect) 

The QB derived depths using the ratio and turbidity algorithms in comparison 
with chart depths are shown in Figs 17 and 18: 


Plymouth Sound Calculated Depths (Ratio Algorithm) 



Figure 17. Ratio algorithm and chart depths of short Plymouth Sound transect. 


Plymouth Sound Calculated Depths (Turbidity Algorithm) 



Figure 18. Turbidity algorithm and chart depths of short Plymouth Sound transect. 
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c. 


ACCURACY COMPARISONS 


1. Looe Key 

The results of derived depths at Looe Key from the ratio algorithm produeed an 
average error of 11.7% along the transeet, with a maximum error of 38%. Fig 16 
highlights the case that as the depth of water increases so to does the inaccuracy of the 
derived depths. In the case of the turbidity algorithm, the average error is 16.5%, the 
maximum error is 36.1%. In the case of the ratio algorithm, the greatest discrepancy in 
depths with the LIDAR data occurs over the darker coral fingers. The turbidity algorithm 
shows the greatest inaccuracy in water deeper than 7m. The algorithm fails to predict the 
increasing depth shown on LIDAR. Another source of inaccuracy in both methods is 
caused by waves and ripples on the sea surface. The higher variability, or noise, in the 
depth profiles in Figs 16 and 17 is, in part, the result of varying reflection of light caused 
by these ripples. In general the ratio algorithm (including the discrepancy over the 
fingers) is more accurate than the turbidity algorithm. 

The QB residual accuracies using the ratio and turbidity algorithms in comparison 
with LIDAR bathymetry are shown in Figs 19 and 20; 


Looe Key Residual Accuracy (Ratio Algorithm) 



Figure 19. Ratio algorithm/LIDAR residual error of Looe Key transect. 
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Looe Key Residual Accuracy (Turbid Algorithm) 



Figure 20. Turbidity algorithm/LIDAR residual error of Looe Key transect. 

2. Plymouth Sound (Short Transect) 

The results from Plymouth Sound show greater variability and noise for both 
methods due to the effects of sea surface variability and higher levels of particulate 
matter in the water column. The amended ratio algorithm average error is 13.1%, with a 
maximum of 28% and a minimum of 0%. The turbidity algorithm average error is 
11.6%, with a maximum of 13% and a minimum of 0.2%. The ratio algorithm does not 
perform as well as the turbidity algorithm in the shallow water, between 3-5m depth. 
However, the turbidity algorithm underestimates depths in excess of 5m. 

The QB residual errors using the ratio and turbidity algorithms in comparison 
with chart depths are shown in Figs 21 and 22: 
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Plymouth Sound Accuracy (Ratio Algorithm) 



Figure 21. Amended ratio algorithm/chart depth residual error of short Plymouth transect. 


Plymouth Sound Accuracy (Turbidity Aigorithm) 



Figure 22. Turbidity algorithm/chart depth residual error of short Plymouth transect. 
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3. Plymouth Sound (Long Transect) 

The effect of increased sedimentary water in Plymouth Sound away from the 
shore is shown in the longer transect in Fig 23. The chart data show steadily increasing 
depth as does the turbidity algorithm to a distance of approximately 350m. At this 
distance a band of more sedimentary water is visible in the QB image (see Fig 14) and 
the resultant calculated depths plateau to give a set of underestimated values of the 
charted depth. 


Plymouth Sound Calculated Depths (Turbidity Algorithm) 



Figure 23. Turbidity algorithm and chart depths of long Plymouth Sound transect. 

4. Sources of Error 

In summary there are several factors that affect accuracy of the QB remotely 
sensed bathymetry calculations. 

a. Effect of Bottom Type 

The coral lingering at Looe Key from approximately 300m to 350m gives 
a deeper set of values in the ratio algorithm rather than shallower as depicted by LIDAR. 
The darker colour of the coral lingers is miss-read by the algorithm as deeper water 
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giving a 38% error in that depth of water. The effeets of bottom type and colour do affect 
the ratio algorithm significantly, and the postulation that the algorithm naturally corrects 
for bottom type is not supported at Looe Key. The differing bottom reflection 
propensities of the QB wavebands alter the integrity of the calculated values. 

b. Effect of Light Atten nation 

The ability of light to penetrate the water column is a limitation for both 
algorithms. The ratio algorithm returns accurate depths to approximately 10m in the case 
of Looe Key, whereas the turbidity algorithm is effective to 7m. In Plymouth Sound, the 
attenuation of light is reduced because the amended ratio algorithm uses the red 
waveband. The absorption of light by water is greater at the red waveband than at the 
blue and green (see Fig 10), and therefore, the effective depth penetration is less. The 
turbidity algorithm bottoms out at approximately 7m as a result. The amended ratio and 
turbidity algorithms reduce noise in the upper water layers because it is scattered less 
than the blue or green wavebands. However, it has a reduced operating depth due to the 
fact that the red waveband is absorbed at a higher rate by the water column than blue or 
green wavebands. 

c. Effect of Surface Waves 

Both methods are affected by the water surface and the noise created by 
variable reflection of light from waves and ripples. The sea surface correction algorithm 
goes a long way to reduce these effects, but does not eliminate them. Figure 24 shows 
the variability of reflection from the sea surface caused by surface ripples, the right hand 
image is taken with the NIR band and gives radiance values for the top most layer of the 
water column. The NIR image is compared with the true colour image, of the same area 
on the left, to highlight how to a lesser effect the visible wavebands are influenced by 
variable surface reflection. Figure 25 gives the Digital Number (DN) values received at 
the satellite for the NIR image and shows the amount of noise caused by sea surface 
ripples. It is assumed that a similar pattern at lower amplitude is present in the red, green 
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Figure 24. True colour and NIR images from the same image and area of Looe Key highlighting 

surface ripple effects on radiance values received at the satellite. 


NIR Raw Values 



Distance (m) 


Figure 25. NIR DN values for the Looe Key transect highlighting the variability of the received 

signal caused by surface ripple effects. 
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and blue wavebands that ereate a sinusoidal pattern in the two algorithms when 
calculating water depth. 

d. Effect of Non-Homogenous Water 

Remote sensing methods rely on the target area having a homogenous 
body of water, both horizontally and vertically. In the case of Plymouth Sound a band of 
run-off crosses the transect with greater amounts of silt and particulate matter. The 
algorithm underestimates depths in the run-off because the increased sediment gives the 
water column a lighter colour, which is miss-read to be shallower than the charted depths. 

D. DEPTH IMAGE OUTPUT 

The ENVI output images for both areas are shown in Figs 26 and 27. These 
bathymetry graphics agree with the chart data for each area and provide a greater degree 
of detail than a standard nautical chart. Small bathymetric changes, due to shifting sands 
and silts for example, can be readily determined through the use of satellite imagery and 
exported in a variety of formats. 



Figure 26. QB derived bathymetry for Fooe Key using the ratio algorithm. 
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Figure 27. QB derived bathymetry for Plymouth Sound using the turbidity algorithm. 
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VI. CONCLUSION AND RECOMMENDATIONS 


A. CONCLUSIONS 

This study uses multi-spectral QuickBird imagery with a nominal resolution of 
2.6m to derive seabed bathymetry at two separate locations, Looe Key (USA) and 
Plymouth Sound (UK). A series of corrections for the atmosphere, the sea surface and 
the water column are employed to convert the blue, green and red wavebands to 
reflectance values. The NPS AOD algorithm, a sea surface glint correction algorithm and 
HYDROLIGHT modeling, is used to determine the corrections required, and the resultant 
reflectance values are incorporated into a ratio (clear water) and turbidity (estuarine 
water) algorithm to derive water depth. The Looe Key results are compared to LIDAR 
bathymetry obtained in 2004 and the Plymouth Sound results compared to hydrographic 
survey data collected between 1963-72. 

The results show that the standard ratio and turbidity algorithms provide an 
accurate representation of seabed bathymetry in the clear water of Looe Key, with 
average errors of 11.7% and 16.5% respectively. In the turbid water of Plymouth Sound, 
the error of the turbidity algorithm averages 11.6%. Although the standard ratio 
algorithm was found not to return valid results at Plymouth, an amended version returned 
results with an average error of 13%. 

The accuracy of QuickBird derived bathymetry does not meet the requirements of 
the International Hydrographic Office’s minimum standards for depth uncertainties for 
inshore waters. However, the low average error of these algorithms, in clear and turbid 
waters, suggests that these methods offer significant potential for determining shallow 
water bathymetry in the littoral zone. 

B. RECOMMENDATIONS 

The performance of both algorithms is dependant on four main environmental 
factors. To improve the accuracy of satellite-derived bathymetry, the following 
recommendations are proposed; 
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1. Resolve Variable Bottom Type 

Variable bottom types in an image should be masked and analysed individually to 
avoid spikes like that over the coral lingers at Looe Key. The target image should be 
separated into areas of similar bottom type and reflection coefficients derived for each 
area to be incorporated into a depth algorithm. 

2. Resolve Light Attenuation 

The extent to which light attenuates in water restricts the depth to which an 
algorithm operates. In turbid water, the QB red and green wavebands are more prone to 
absorption, and therefore, are limited in the depth to which they can retrieve accurate 
depth estimates. Similar depth analyses in turbid water with different sensors that have 
different red and green waveband sensors may provide greater depth penetration than 
those of QB. 


3. Resolve Surface Effects 

An analysis of surface waves and ripples could be used to determine a sinusoidal 
wave pattern in the derived depth profile. If this was subtracted from the depth profile, or 
a dedicated sampling procedure developed, the noise could be corrected for. 

4. Resolve Non-Homogeneous Water 

An analysis to determine the horizontal variability of matter in the water column, 
such as the chlorophyll method developed in this study, would give more accurate water 
coefficient results when modeled in software such as HYDROLIGHT. The results when 
incorporated into HYDROLIGHT as individual pixel values would give attenuation 
coefficient values for each pixel. Such values will enable an algorithm to analyse an area 
of non-homogeneous water more accurately. 
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